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El problema del análisis de la flexión simple en piezas rectas de sección constante con ciertas 
hipótesis complementarias y otras condiciones de forma y solicitación nos conduce a un modelo regido por 
una ecuación diferencial lineal de coeficientes constantes con condiciones de contorno que puede ser 
reducido a un problema unidimensional definido en un intervalo. Este problema es, como se sabe, 
fundamental en el campo de la edificación por el predominio de estas piezas en arquitectura e ingeniería. 
La presencia de cargas y momentos puntuales en las solicitaciones de vigas es imprescindible para 
completar un modelo operativo ya que permite ampliar substancialmente el campo del análisis a vigas de 
varios tramos y otros problemas de interés en el contexto estructural. Pero esta presencia nos impide 
efectuar holgadamente integraciones "convencionales" con lo que se necesita acudir a otros procedimientos 
matemáticos entre los que se encuentran los métodos numéricos. 

El método de integración por diferencias finitas es un método numérico universal y 
extraordinariamente eficaz que resulta una solución idónea por su simplicidad y adaptabilidad a cualquier 
tipo de planteamientos en este ámbito. El sistema de ecuaciones al que conduce la resolución de los 
problemas puede hoy día ser resuelto con gran facilidad por la mayor parte de calculadoras ordinarias y 
programas informáticos. 

Después de una exposición del procedimiento con carácter general descrito sobre los casos más 
simples para ecuaciones de 2” y 4” orden, se ha procedido a configurar unas fórmulas de uso sistemático 
totalmente adaptadas a las condiciones del problema de la flexión simple en lo que hemos denominado 
"métodos en dos etapas” partiendo de la caracterización de las condiciones de conjugación de la unión 
rígida, el tratamiento de las secciones extremas y, en su caso, de la unión articulada. Con esta descripción 
pormenorizada puede ya procederse reiteradamente sobre las piezas de uno o varios tramos sean cuales 
sean sus vinculaciones de extremo. Se ha dado la mayor importancia a los ejercicios prácticos que son los 
que permiten comprender la metodología operativa y obtener los resultados buscados en los distintos casos 
comparando sus soluciones con las obtenidas por otros métodos. 

Finalmente somos conscientes de que en un trabajo de estas características es inevitable la aparición de 
erratas o de falta de precisión que iremos corrigiendo en revisiones posteriores. 


Pedro Galán del Sastre 
Ramón J. Zoido 
Mayo 2014 
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2 Solicitaciones segmentariamente 
continuas con saltos finitos 


2.1. La flexión simple en vigas rectas de rigidez 
constante 


Es sabido que el análisis de la flexión simple de barras rectas de rigidez constante es 
fundamental en el campo de la edificación por el predominio de estas piezas en arquitectura 
y, en general, también en ingeniería. Aquí trataremos un conjunto de problemas que pueden 
ser resueltos de una manera muy simple y muy eficaz con el auxilio de métodos numéricos en 
diferencias. Soslayando conceptos tales como lo ¿sostático o lo hiperestático —conceptos que, 
por otra parte, pueden deducirse aquí del conjunto final de ecuaciones resultante—, se procede 
de manera sistemática sobre todo el rango de piezas sean cuales sean sus vinculaciones de 
extremo aunque aquí plantearemos los casos estándar más habituales de extremos libres, 
simplemente apoyados o perfectamente empotrados y sus posibles combinaciones siempre 
que determinen modelos mecánicamente estables. Con el mismo tratamiento, el método se 
extiende a las vigas continuas incluyendo los apoyos interiores rígidos— y las compuestas 
incluyendo las articulaciones internas-. 

Las vigas que vamos a considerar aquí son piezas horizontales de directriz recta —eje OX 
longitudinal- y sección transversal constante con un eje de simetría —paralelo al OY-. Si la 
longitud de esta pieza l resulta ser una magnitud mucho mayor que las dimensiones de la 
sección, podemos considerar que la flexión queda definida por la deformación de la sección 
longitudinal media de la viga, en el plano de la directriz OX Y —plano de simetría de la 
viga— cargada verticalmente en este plano con una ley de cargas cuya variación viene descri- 
ta por la función w (1). Admitiremos además, por una parte, la hipótesis de existencia de 
fibra neutra —fibra longitudinal sin extensión ni contracción— y, por otra, consideraremos que 
las secciones transversales inicialmente planas y perpendiculares a la línea neutra lo siguen 
siendo después de la deformación (Bernouilli-Navier). Esta última hipótesis desprecia en la 
flexión el alabeo de la sección producida por el esfuerzo de corte. Básicamente con estas 
hipótesis y estableciendo el origen en el extremo izquierdo de la pieza, el problema de la 
deformación puede ser reducido a un problema unidimensional definido en el intervalo |[0, /] 
(figura 2.1). Los esfuerzos se limitan a los momentos flectores y a los esfuerzos de corte o 
cortantes en tanto que las deformaciones se reducen al desplazamiento vertical — flechas— y 
al giro de cada sección respecto de su posición inicial —giros—. Es muy significativo el modelo 
de secciones transversales rectangulares. Siendo el “ancho” de cada sección transversal una 
magnitud pequeña, podemos considerar que todas las secciones longitudinales se comportan 
idénticamente en la deformación. La deformada de la línea neutra y (1) que, en estas condi- 
ciones, domina el problema de la deformación, se denomina curva elástica. Las secciones 
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2 Solicitaciones segmentariamente continuas con saltos finitos 
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Figura 2.1: Reducción de dimensión del problema de la deformación. 


transversales permanecen perpendiculares a esta línea y su giro puede ser determinado por 
el ángulo que forma la recta tangente a esta curva con la horizontal que es la posición de la 
fibra neutra antes de la deformación. Esta curva deformada viene regida por la ecuación 


y" = 5 (1 de y?) 


|w 


(2.1) 


donde m = m (1) es el momento flector que solicita la pieza y que puede ser descrito física- 
mente en un punto de la barra como el momento —de primer orden— de las fuerzas situadas a la 
izquierda del punto incluyendo las de reacción de su extremo izquierdo 
—reacciones y momentos— si las hay. A su vez B factor de rigidez a flexión de la viga es 
igual al producto El donde E es el módulo de elasticidad de Young que depende del tipo 
de material y que supondremos constante a lo largo de la pieza e [ = I, es el momento 
de inercia de la sección respecto de un eje normal al plano de la directriz por el centro de 
gravedad de la misma. La ecuación (2.1) es válida en el caso más general de vigas de sección 
variable B = B(x) a lo largo del eje longitudinal y manifiesta simplemente la relación de 
proporcionalidad local entre la curvatura de la curva elástica y los momentos flectores siendo 
justamente el factor de proporcionalidad el valor B de la rigidez a flexión. 

El dimensionado de la pieza y las limitaciones en la ley de cargas en cada caso garantizan 
que la función y (1), solución de esta ecuación de campo más las condiciones de extremo, 
es continua en el intervalo (0,/) y es al menos de clase uno si no existen en dicho intervalo 
articulaciones internas que conviertan a estas piezas en vigas compuestas. Una discontinuidad 
“natural” en la primera derivada de la viga supondría una fisura semejante a una rótula y, 
en cierto modo, su alejamiento del rango elástico de respuesta o incluso su alejamiento del 
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2.1 La flexión simple en vigas rectas de rigidez constante 


comportamiento de viga pudiendo transformarse en un mecanismo. 

Si consideraremos que en la deformación los giros resultantes de las secciones son pequeños 
podemos trabajar con un modelo lineal aproximado para la ecuación de campo ya que sien- 
do y” pequeño, el valor y? puede ser depreciado. Esta hipótesis es plausible toda vez que 
giros importantes de las secciones serían inviables con los materiales utilizados usualmente y 
mecánicamente constituirían otra respuesta diferente de la flexión simple por lo que deberían 
ser abordados con otros procedimientos. La consideración 1+y” = 1 es, por otra parte, con- 
sistente con la propia hipótesis de la existencia de la fibra neutra —fibra sin alargamiento 
ni contracción— ya que este factor es justamente el cuadrado del integrando que determina 
la longitud de la curva elástica. Finalmente esta hipótesis es congruente con los resultados 
experimentales obtenidos para las piezas que habitualmente se utilizan en construcción. 

Con esta hipótesis adicional la ecuación de campo que rige la flexión simple de las vigas 
rectas en el modelo lineal correspondiente resulta ser simplemente: 


y" (1) = a 


Si consideramos, además, vigas de inercia constante, [, = cte, y con el módulo de elasti- 
cidad de Young constante, E = cte, el modelo lineal satisface la ecuación 


con B = El = cte, y constituirá la ecuación de campo de las vigas de sección —o rigidez— 
constante. 
En estas condiciones, la derivación de esta ecuación nos lleva a 


m 
donde t (x) es el esfuerzo de corte o cortante, y una nueva derivación 
] w (2) 
10] = 2.2 
y)= 22) 


donde w (1) es la ley de cargas que solicita la pieza a lo largo del intervalo |[0, /]. 

La ecuación (2.2) es la ecuación diferencial de la curva que representa los descensos de 
cada una de las secciones transversales de la pieza en las condiciones señaladas y resulta 
ser una ecuación de cuarto orden lineal cuya simplicidad dependerá de la función w (x) del 
segundo miembro. Si podemos integrar directamente esta ecuación para obtener y (1), flecha 
o desplazamiento de cada sección, entonces podemos completar la deformación de la viga 
obteniendo por derivación en primer lugar 


y (u) =tg0 (1) = 0 (2), 


giro de cada una de las secciones transversales en cada punto de la viga, mientras que las 
siguientes derivadas nos proporcionarán los esfuerzos de flexión y corte: 


By le] = miz), momento flector, 


By (2) = 10), esfuerzo cortante. 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


2.2. Las leyes de cargas w (1) 


Consideraremos aquí tres tipos de acciones verticales que representarán un amplio rango 
de solicitaciones generales: 


a. Una distribución de cargas verticales definida por una función segmentariamente con- 
tinua con un número finito de discontinuidades de primera especie —saltos finitos— en 
el intervalo. 


b. Un número finito de cargas concentradas en distintos puntos del intervalo o cargas 
“puntuales”. 


c. Un número finito de momentos de flexión aislados definidos en distintos puntos del 
intervalo. 


Cuando trabajamos solamente con el primer conjunto “a” de cargas, podremos abordar direc- 
tamente la ecuación de cuarto orden y realizar hasta donde nos sea posible el cálculo discreto 
de las funciones de deformación y esfuerzo en lo que denominaremos método de una etapa en 
el cual básicamente obtenemos por métodos en diferencias las flechas de la viga aproximando 
la integración de la derivada de orden cuatro y de ahí reconstruimos por derivación el resto 
de las funciones —giros y funciones de esfuerzo—. La consideración del segundo y tercer grupo 
“b” y “c”, nos lleva preferentemente a articular los métodos en dos etapas con la introducción 
de las cargas concentradas y los momentos aislados a través de las condiciones de conju- 
gación de la unión rígida y el replanteamiento de las condiciones de los extremos estándar. 
La presencia de estos dos grupos de cargas “b” y “c” nos permitirá, además, extender estos 
métodos a las vigas de varios tramos. 


2.3. Condiciones de contorno. Aproximación de las 
condiciones de contorno en ausencia de cargas 
concentradas y momentos aislados 


La descripción pormenorizada de la caracterización de los vínculos extremos en las vigas 
cae dentro del análisis del grado de empotramiento de las estructuras pero aquí considera- 
remos un conjunto de condiciones de frontera estándar, idealizadas desde el punto de vista 
estructural, que pueden caracterizarse de forma muy simple pero que permiten completar 
modelos de magníficos resultados con respecto al comportamiento real. En ausencia de cargas 
puntuales y momentos aislados en los extremos serán las siguientes: 


a. Extremo libre: caracterizado por la anulación de las derivada segunda —momento flector 
nulo, giro libre- y tercera cortante nulo-. 


b. Extremo simplemente apoyado: caracterizado por la anulación de la flecha —apoyo 
rígido, no hay descenso— y de la derivada segunda —momento flector nulo, giro libre sin 
restricciones angulares-. 
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2.3 Condiciones de contorno. Aproximación 


tipo de viga extremo izquierdo extremo derecho 


Y. 


ménsula empotrada en su 
extremo izquierdo 


a 


ménsula empotrada en su 
extremo derecho 


¿¿g Ñq----——— A 


viga doblemente apoyada 


Ha 


viga empotrada-apoyada 


| 


viga apoyada-empotrada 


J l 


viga doblemente empotrada 


Tabla 2.1: Leyes de cargas segmentariamente continuas sin acciones aisladas en los extremos. 


c. Extremo perfectamente empotrado: caracterizado por la anulación de la flecha 
-empotramiento rígido, no hay descenso- y de la derivada primera -giro nulo de la 
sección de empotramiento, la recta tangente a la curva elástica es horizontal-. 


Las configuraciones resultantes de combinar estas tres condiciones de extremo en un solo 
tramo de manera que sean mecánicamente estables y que proporcionen reacciones resistentes 
de viga, conducen a los seis modelos de sustentación que se muestran en la tabla 2.1. A la 
derecha de cada uno de los seis modelos se representan sus características en términos de 
los valores de la curva elástica o de sus derivadas. Estas condiciones son válidas en ausencia 
de cargas concentradas y momentos aislados en los extremos. Cuando, además, se conside- 
ren solicitaciones puntuales en estos extremos —idealmente representadas en los esquemas 


Figura 2.2 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


vínculos condiciones de extremo 
empotramiento perfecto yo =0 Yo =0 
+- HA 
1 12 9>Y1_7 aan 
| 2h sil 1 
apoyo rígido Yo =0 Yo =0 
“2 
“1 12 n= Yo tYr_y Ñ 
p2 7. Y-1 — Yi 
—= 2Y0 + Y- 
extremo libre el E 8 0 Y-1 = 2Y0 — Yi 
+ 4 44 
-2-10 12 Ya — 2y1 + 2y-1 — y- 
- an - 2=0 | y-2= Yo — 4Y + Ya 


Tabla 2.2: Solicitaciones segmentariamente continuas sin acciones aisladas en los extremos. 


de la figura 2.2, o sea: cargas concentradas en bordes libres y/o momentos aislados sobre 
bordes libres o apoyos rígidos—, las condiciones de la tabla deberán ser convenientemente 
modificadas. La introducción de este tipo de solicitaciones se verá en el $ 3. 

En el caso en que las solicitaciones no incluyan cargas concentradas ni momentos aislados 
entonces teniendo en cuenta las aproximaciones O (h?) que venimos considerando, las condi- 
ciones de extremo estándar de la tabla 2.1 que vienen descritas en términos de derivadas 
permiten asignar a los puntos “ficticios” ciertos valores de la función buscada respecto de los 
que toma la función en los nodos del interior del intervalo más próximos al contorno. 

Tomando el extremo izquierdo 1y = 0 y dos nodos interiores 21 y 12 Como referencia, se 
muestran en la tabla 2.2 las condiciones de extremo y los valores que se deben de asignar a los 
puntos ficticios x_, y, en su caso, 1_2 según las condiciones del vínculo respecto a los valores 
de los puntos interiores próximos al punto frontera xp. Estas condiciones se completan en la 
tabla con Yo = 0, descenso nulo en apoyos rígidos y empotramientos perfectos. 


2.4. Métodos de una etapa. Ejemplos 


En ausencia de cargas concentradas y momentos aislados se puede abordar directamente 
la ecuación de cuarto orden (2.2). Aproximaremos la integración de la ecuación (2.2) más las 
condiciones en los extremos teniendo en cuenta el factor de escala B = cte que debe reducir 
(por división) los resultados de esta integración directa para los valores de las flechas y giros 
en tanto se conservan los resultados sin reducción de escala en las funciones de esfuerzos. El 
procedimiento en diferencias para este tipo de ecuaciones ha sido descrito ya en $ 1. 

Los siguientes ejemplos ilustran cada uno de los seis casos de la tabla 2.1 sometidos a 
leyes de cargas representativas más o menos asimilables a tipos de cargas reales que pueden 
solicitar las piezas. En todos los casos haremos coincidir las posibles discontinuidades sobre 
los nodos de la red y, en su caso, cuando haya que evaluar la ley de cargas en los extremos será 
a la derecha o la izquierda del nodo extremo siempre del lado interior del intervalo. En cada 


36 


2.4 Métodos de una etapa. Ejemplos 


uno de ellos, obtendremos los diferentes conjuntos de soluciones discretas sobre tres redes 
progresivamente más finas de pasos h =1m,h=0.5 m y h= 0.25 m y los compararemos 
con los resultados exactos calculados por métodos continuos. Obtenida la proyección discreta 
de los descensos o flechas de las secciones para h = 0.25 m, reconstruiremos a partir de estos 
valores, por derivación, sus giros y los valores de las funciones de esfuerzo completando así 
el análisis de cada una de las vigas. 


2.4.1. Ménsula empotrada en su extremo izquierdo. 


Se considera la viga de 4 m en ménsula —empotrada-libre- con sobrecargas uniformes 


definidas en la forma: 
js 2, si0<zxw<2, 
E 


tal y como muestra la figura 2.3(a). 

Estableceremos en primer lugar una red discreta de paso h = 1 m y enumeraremos los 
nodos en la forma que se representa en la figura 2.3(b). Aplicando la ecuación de campo, con 
h = 1, sobre cada nodo interior (o extremo) y teniendo en cuenta los valores de los nodos 


(b) 
2 2 2 
3 4 
A AAA 3) ==2222218y=4+= Aer + 
-1 ÉS 1 2 5 6 
0 1 1 1 
(c) 
2 2. 02. 2 2 
506 7 8 
O E OS O AS SS A a de. ---- 
ES ES 1 2.34 o 10 
1.1. 1 1.1 


(d) 
1121222222 


9 1011 12 13 14 15 16 
+ ++ 
1 1234567? 7 18 


111111111 


Figura 2.3 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


exteriores asignados a través de las condiciones frontera, el sistema de ecuaciones resultantes 
será: 


Y-1—4AYo +6Y1 —4Ya+Yg = 2, con y-1 =Y1, Yo =0 

Yo — 4y1 + 6y2 — 4ya + Ya = 0.5, con yo =0 

y — 4ya+6y3 — Aya+ ys = —l, con Ys = 2Ya — Y3 

Ya — 4y3 + Oya — 4Ys + Ys = —l, con Ys = 2Ya — Y3, Yo = Aya — Ay3 + Ya 


que matricialmente se puede escribir como 


7 -4 1 0 Yi 2 
4 6 -4 1 Y | |05 
1-4 5 -2 y | | -1 
0 2 -4 2 Ya =1 
y cuya solución es 

Yo 0.0 

Yi —1.0 

Ya = —5.0 


Como ya vimos en $ 1, podemos escribir el sistema de forma matricial para un valor de h 
genérico, lo que nos permitirá resolver para cualquier valor de h: 


Fo od 1 Yi w(x1) 
2 € = e 1 Ya w (22) 
1-46 -4 1 Ya w (23) 
5 1. -46 4 1 Yi = w (x;) 
E o | UYN-2 w(2y-2) 
LA =L UYN-1 w(2y-1) 
a UN w (ty) 


De esta forma, para obtener una mejor aproximación, consideraremos ahora una red de 
paso h = 0.5 sobre la notación de la figura 2.3(c), y otra de paso h = 0.25 figura 2.3(d)-—. Los 
resultados obtenidos para los distintos valores de h pueden ser comparados entre si y con los 
valores exactos obtenidos por métodos continuos en la tabla 2.3. Las poligonales asociadas 
a los valores discretos pueden compararse con los exactos en la figura 2.4. 

Eligiendo los resultados de mejor aproximación, los obtenidos sobre la red Ry-0.25, pueden 
ser calculados ahora los momentos flectores en los puntos de la red con la fórmula: 

— Yi+1— 2Yi + Yi 


m¡= 0.952 , con ¿=0,1,2,...,16, 


donde y-¡ e yi7 se pueden calcular mediante las fórmulas obtenidas en la tabla 2.2 y las 
análogas para los puntos ficticios a la derecha del intervalo. 
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2.4 Métodos de una etapa. Ejemplos 


x discretos continuos 
h= h=0.5|h=0.25 

0.00 0.0000 0.0000 0.0000 0.0000 
0.25 -0.0625 -0.0674 
0.50 -0.2500 | -0.2773 -0.2865 
0.75 -0.6641 -0.6768 
1.00 -1.0000 | -1.1875 -1.2344 -1.2500 
1.25 -1.9922 -2.0101 
1.50 -2.87150 | -2.9336 -2.9531 
1.75 -4.0469 -4.0674 
2.00 -5.0000 | -5.2500 | -5.3125 -9.3339 
2.25 -6.7031 -6.7241 
2.50 -8.1250 | -8.1895 -8.2109 
2.15 -9.7461 -9.7684 
3.00 || -11.0000 | -11.2812 | -11.3516 -11.3750 
3.29 -12.9883 -13.0132 
3.90 -14.5625 | -14.6426 -14.6693 
3.19 -16.3047 -16.3339 
4.00 | -17.5000 | -17.8750  -17.9687 -18.0000 


Tabla 2.3: Resultados obtenidos, By, para distintos valores de h y solución exacta. 


Para completar el análisis de las funciones de deformación y esfuerzo sobre este conjunto 
de resultados correspondiente a la red h = 0.25, podemos obtener los giros de cada una de 
las secciones en los nodos de la red sin más que realizar las operaciones: 


0 rn Yi41 — Yi-1 


¡== Y, = E j con+=0,1.2..25 10, 


sin olvidar que hay que dividir los valores obtenidos por el factor B si queremos disponer de 
estos giros en su verdadera escala. 
Finalmente puede construirse en forma discreta, hasta donde es posible, la ley de cortantes 


5 


-10 


-15 


1 2 3 


Figura 2.4: Solución analítica junto con la solución numérica para h=1,h=0.5 y h = 0.25. 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


Figura 2.5: De arriba a abajo: ley de cargas By?”, cortante By”, flectores By”, giros By” y flechas By. 


a través de las fórmulas: 
mm Yi+2 — 2Yi41 + 2Yi-1 — Yi-2 
E 2 - 0.253 
lo que nos proporciona directamente los esfuerzos de corte en cada uno de los nodos señalados. 
En el apartado $ 3.5 completaremos estos cálculos de manera que incluyan estos esfuerzos 
sobre los vínculos. 

La reconstrucción de las gráficas correspondientes a los valores discretos de flechas, giros, 
momentos y cortantes se muestran finalmente en la figura 2.5 bajo la correspondiente ley de 
cargas. 


ti=Yy con ¿=1,2,...16, 


2.4.2. Ménsula empotrada en su extremo derecho 


Consideramos una viga de 4 m en ménsula, perfectamente empotrada en su extremo 
derecho figura 2.6(a)- cuya ley de cargas puede ser descrita en la forma: 


ES si0<x<2, 


US ateancd 


Establecemos en primer lugar una red discreta de paso h = 1 m y enumeramos los nodos 
en la forma que se representa en la figura 2.6(b). El sistema que obtenemos será: 


Y-2 — Ay-1+6y0 —4y1+Y2a = 0, con Y-2 = 4Yo — 4Y1 + Ya, Y-1 = 2Yo — Y1 
Y-1—AYo+6y1 — Ay + Ya = 4, con Yy-1 = 2Yo — Y1 

Yo — Ay1 + 6y2 — Aya +Yya = 5, con ya =0 

y — Aya + Oya — MYya+Ys = 2, con Ya =0, Ys = Y3 


40 


2.4 Métodos de una etapa. Ejemplos 


++- 
-2-1012345678 9101112131415 417 


Figura 2.6 


que matricialmente se puede expresar como: 


2 -—4 2 0 YO 2 
2 5 -4 1 Y | 1.05 
1-4 6 -4 Y | | -1 
0 1 -4 7 Y3 =1 
y cuya solución es 

Yo 95.0 

Yi 66.0 

y |=| 37.0 

Y3 12.0 

Ya 0.0 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


x discretos continuos 
h=1|h=05|h=0.25 

0.00 | 95.0000 | 98.0000 | 98.7031 98.9333 
0.25 90.7500 90.9334 
0.50 82.3790 | 82.7969 82.9344 
0.75 74.8477 14.9412 
1.00 | 66.0000 | 66.7500 | 66.9141 66.9667 
1.25 99.0195 99.0351 
1.50 91.2500 | 51.2031 91.1865 
1.75 43.5234 43.4804 
2.00 | 37.0000 | 36.2500 | 36.0625 36.0000 
2.25 28.9297 28.8945 
2.50 22.5000 | 22.2539 22.1719 
2.15 16.1719 16.0889 
3.00 | 12.0000 | 11.0625 | 10.8281 10.7500 
3.20 6.3790 6.3076 
3.90 3.1250 2.9727 2.9219 
3.19 0.7891 0.7607 
4.00 0.0000 | 0.0000 0.0000 0.0000 


Tabla 2.4: Resultados obtenidos, By, para distintos valores de h y solución exacta. 


Igual que en el ejemplo anterior, escribimos el sistema para un valor h genérico de forma 


matricial: 


2 


—4 
4 6 


4 1 
4 6 -4 1 
1-46 -4 1 
4 6 -4 
1. -4 7 


Yo 
Yi 
Ya 


Yi = 


UN-3 
YN-2 
YN-1 


y resolvemos para una retícula de paso h = 0.5 figura 2.6(c)- y h = 0.25 figura 2.6(d)-. 
Los resultados —que deben dividirse por el factor B para obtener los correspondientes 
descensos— para los distintos valores de h pueden ser comparados entre si con los valores 
exactos obtenidos por métodos continuos en la tabla 2.4 y así mismo podemos comparar las 
gráficas obtenidas para las tres discretizaciones con los resultados exactos en la figura 2.7. 


Eligiendo los resultados obtenidos para la proyección sobre la red Ry-=0.25, podemos hallar 
los momentos flectores en los puntos de la red con los cálculos complementarios: 
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Yi 2Ui PU 


2 


(1.208 


») 


con +=0,1,2,...,16; 


2.4 Métodos de una etapa. Ejemplos 


—S— exacta 
——- h=0.5 


Figura 2.7: Solución analítica junto con la solución numérica para h=1,h=0.5 y h = 0.25. 


los giros de cada una de las secciones: 
rn Yi+1— Yi-1 
[== RU con¿+=0,1,2,...,16; 
da 0.5 
y las cortantes en cada uno de los nodos a través de las fórmulas: 
m Visa — 2Yig1 + 2Yi-1 — Yi-2 
es 20.25 | 
La reconstrucción de las gráficas correspondientes a los valores discretos de flechas, giros, 
momentos y cortantes se muestran finalmente en la figura 2.8 bajo la correspondiente ley de 


con¿=0,1,2,...,15. 


Cargas. 


Figura 2.8: De arriba a abajo: ley de cargas By", cortante By”, flectores By”, giros By” y flechas By. 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


2.4.3. Viga doblemente apoyada 


Aquí consideraremos la viga de 6 m con dos apoyos rígidos en sus extremos (figura 2.9(a)) 
solicitada por una ley de cargas que puede venir representada por la función: 


—=2+4, si0<zx<2, 
wla)=x 1, 12<6 4, 
=x+ 7, si4<x<6. 


Podemos iniciar nuestro proceso de aproximación estableciendo una red de paso h = 1 m. 


7 
710203040506 0730910111301 015 16 17 18 19 20 21 2 2225 


Figura 2.9 
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2.4 Métodos de una etapa. Ejemplos 


de discretos continuos 
h=1|h=05|h=0.25 

0.00 0.0000 | 0.0000 0.0000 0.0000 
0.25 4.0677 4.0605 
0.50 8.0747 8.0394 8.0222 
0.75 11.8179 11.7997 
1.00 | 15.6667 | 15.4097 | 15.3435 15.3213 
1.25 18.5534 18.5279 
1.50 21.4844 | 21.4004 21.3721 
1.75 23.8479 23.8175 
2.00 | 26.3333 | 29.96593 | 25.8694 25.8370 
2.25 27.4468 27.4127 
2.50 28.6753 | 28.5681 28.5322 
2.15 29.2253 29.1877 
3.00 | 30.0000 | 29.5312 | 29.4141 29.3750 
3.25 29.1341 29.0939 
3.90 28.5122 | 28.3889 28.3480 
3.19 27.1860 27.1450 
4.00 | 26.1667 | 25.6597 | 25.5369 25.4963 
4.25 23.4568 23.4177 
4.50 21.0781 | 20.9688 20.9326 
4.75 18.1067 18.0745 
5.00 | 15.3333 | 14.9965 | 14.9143 14.8870 
5.25 11.4438 11.4225 
9.90 7.8003 7.7996 7.7408 
9.19 3.9167 3.9091 
6.00 0.0000 | 0.0000 0.0000 0.0000 


Tabla 2.5: Resultados obtenidos, By, para distintos valores de h y solución exacta. 


Con la numeración de los nodos de la figura 2.9(b), obtenemos el sistema de ecuaciones: 


Y-1—4AYo +6Y1 —4Ya+Yg = 3, con y-1 = —Y1, Yo =0 
Yo — 4y1 +6y2 —4ya+Yya = 1.5, con yo =0 

y —4y2+6y3—4Aya+Yys = 1, 

Ya — 4y3 + 6ya — 4Ys + Ys = 2, con Ys =0 

Ya — Aya + 6ys — Aye + yr = 2, con Yé6 = 0, y7 = —Ys 


que matricialmente expresamos como: 


| 0 Y 3 
cd o de 10 Ya 1.5 
LD ld y |=| 1 
0 1 -4 6 1 Ya 2 
0. 0 1 -45 Ys 2 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


Figura 2.10: Solución analítica junto con la solución numérica para h= 1, h=0.5 y h = 0.25. 


y cuya solución es 


Y6 


0.0000 
15.6666 
26.3339 
30.0000 
26.1666 
15.3333 

0.0000 


En este caso, el sistema para un valor h genérico viene dado, en forma matricial, por: 


5 —á4 


4 6 


1 -4 6 


—4 


1 


Yi 
Ya 
Y3 


Yi 


-4 6 -4 1 dra 
> 0 dé YN-2 


4 5 YN-1 


Estableciendo ahora un paso de red h = 0.5 con la notación de la figura 2.9(c) y una red 
discreta de paso h = 0.25 con la representación de la figura 2.9(d), resolvemos nuevamente 


el sistema obtenido. 


Los resultados para los distintos valores de h pueden ser comparados entre si y con los 
valores exactos obtenidos por métodos continuos en la tabla 2.5. También podemos comparar 
las gráficas construidas con las poligonales para los tres valores del paso de red con los valores 


exactos en la figura 2.10. 


Eligiendo los resultados obtenidos en la proyección sobre la red Ry=0.25, pueden ser calcu- 
lados los momentos flectores en sus puntos a través de los cálculos complementarios: 
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Yi 2Yi FB 


(1.258 


cons=0.1.2.34522, 


2.4 Métodos de una etapa. Ejemplos 


Figura 2.11: De arriba a abajo: ley de cargas By?”, cortante By””, flectores By”, giros By” y flechas By. 


y, análogamente, obtener los giros de cada una de las secciones con 


rn Yi+1 — Yi-1 


Y; 05 4 con+=0,1,2,...,24. 


Para obtener los esfuerzos cortantes en esta viga, habrá que aplicar las fórmulas de aproxi- 
mación: 


e Yia — LYa + ia — 2 


; 2.025 : conti= 1,2... 20: 


Estos resultados junto a los momentos, giros de las secciones y flechas para la aproximación 


h = 0.25 pueden verse sobre las gráficas de las figura 2.11 bajo la ley de cargas que solicita 
la pieza. 


2.4.4. Viga empotrada-apoyada 


Se considera aquí la viga empotrada-apoyada de 5 m representada en la figura 2.12(a) y 
sometida a cargas triangulares definidas por la función continua: 


(0) —=2+2, sil0<zx<2, 
w(1)= 
z—2, si2<x<5. 


Estableciendo en primer lugar la red de paso h = 1 m con la notación de la figura 2.12(b), 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


obtenemos el sistema de ecuaciones: 


Figura 2.12 


Y-1 — 4Yo + 6y1 — 4Ya + Y3 


Yo — 
Yi 


Ya 


4y1 + 6ya — 4y3 y 
4y2 + 6y3 — 4Ya - 
4y3z + 6ya — 4ys 


F Ya 
F Ys 


Y6 


con Y-1 =Y1, Yo =Ú 
con Y =0 
con Ys =0 
con Ys = 0, Ys = —Ya 


que matricialmente se puede expresar de la siguiente manera: 
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7 -4 1 0 
4 6 -—4 1 
1-46 -4 
0 1 -—4 5 


N RO 


2.4 Métodos de una etapa. Ejemplos 


x discretos continuos 
h=1|h=05|h=0.25 

0.00 | 0.0000 | 0.0000 0.0000 0.0000 
0.25 0.0919 0.0848 
0.50 0.3996 0.3212 0.3092 
0.75 0.6482 0.6334 
1.00 | 1.2333 | 1.0858 1.0395 1.0232 
1.25 1.4660 1.4494 
1.50 1.9475 1.9029 1.8868 
1.75 2.3282 2.3130 
2.00 | 2.8471 | 2.7603 2.1219 2.7076 
2.25 3.0649 3.0912 
2.50 3.3148 3.3382 3.3247 
2.19 3.9237 3.9104 
3.00 | 3.7412 | 3.6418 3.6053 3.9921 
3.20 3.9697 3.0967 
3.90 3.4433 3.4078 3.3992 
3.19 3.1151 3.1031 
4.00 | 2.8233 2.7239 2.6930 2.6821 
4.25 2.1498 2.1406 
4.50 1.5217 1.5017 1.4948 
4.75 0.7736 0.7697 
5.00 | 0.0000 | 0.0000 0.0000 0.0000 


Tabla 2.6: Resultados obtenidos, 


y cuya solución es: 


El sistema escrito de forma matricial para un paso genérico h queda: 


7 -—4 1 
4 6 -4 1 
1-4 6 -4 


== 1  -—A 


il 


By, para distintos valores de h y solución exacta. 


0.0000 
1.2353 
2.8471 
3.7412 
2.8235 
0.0000 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


—e— 


exacta 
h=1 


exacta 


—e— 


h=0.5 


exacta 
—8— h=0.25 


1 2 


3 


2 


3 


2 


3 4 5 


Figura 2.13: Solución analítica junto con la solución numérica para h= 1, h=0.5 y h = 0.25. 


Así, resolvemos para una red de paso h = 0.5 m, con la notación de la figura 2.12(c) y 
para una red de paso h = 0.25 m con la numeración de nodos de la figura 2.12(d). 


Los resultados para los distintos valores de h pueden ser comparados entre si y con los 
valores exactos obtenidos por otros métodos en la tabla 2.6. Pueden visualizarse estas com- 
paraciones con las gráficas de las poligonales asociadas a estos tres conjuntos de soluciones 
en la figura 2.13. 


Eligiendo los resultados obtenidos para la proyección sobre la red h = 0.25, pueden ser 


Figura 2.14: De arriba a abajo: ley de cargas By?”, cortante By””, flectores By”, giros By” y flechas By. 
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2.4 Métodos de una etapa. Ejemplos 


calculados los momentos flectores con los cálculos complementarios: 


1 = 2Y; i— : 
img = Y A Lo coni=0,1,2,...20. 


Para obtener los giros de las secciones se aplicarán las fórmulas: 


y os Vi+1 7 Yi-1 
y 0.5 : 


Finalmente para obtener los esfuerzos de corte aplicaremos: 


con ¿+=0,1,2,...,20. 


¡+2 — 2Y; 2Yi-1 — Yi 
¡ea te Yi + áYia Y A con?= 120. 10 
2 (1.258 
La ley de cargas, los esfuerzos de corte, los momentos, los giros de las secciones y las flechas 
halladas para la aproximación h = 0.25 se representan gráficamente en la figura 2.14. 


2.4.5. Viga apoyada-empotrada 


Viga de 6 m apoyada-empotrada sometida a una ley de carga triangular como se muestra 
en la figura 2.15(a) que puede venir descrita analíticamente por la función: 


le De si0<zx<3, 
wi = 
—2+6, si3<zx<06. 


Estableciendo como retícula de paso h = 1 representada en la figura 2.15(b), se puede 
obtener el sistema de ecuaciones: 


Y-1— 4Y0+6Y1 —4y9+Ys = 1, con Yy-1 = —Y1, Yo =0 


j 
Yo — Ay + 6y2— Aya+ Ya = 2, con yo =0 
Yi — Aya + Óya = Mya + Ys = 3, 
Ya — 4y3 + Oya — 4Ys + Y6 = 2, con Ya =0 
Ya — 4Ya + 6Yys — Aye + Yy7 = 1, con Yó6 = 0, Y7 = Ys 


que matricialmente podemos expresar como: 


5 -4 1 0 0 Yi 1 
4 6 -4 1 0 Ya 2 
1-4 6 -4 1 y |=]| 3 
0 1 -4 6 -4 Ya 2 
0. 0 1 -4 7 Y5 1 
y cuya solución es: 

Yo 0.0000 

Yi 8.9795 

Ya 14.8767 

y | =| 15.6096 

Ya 11.0959 

Y5 4.2534 

Y6 0.0000 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


78 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 


Figura 2.15 


Si escribimos el sistema de forma matricial para una retícula de paso h: 


5 -4 1 yr w (21) 
4 6 —A 1 Ya w (27) 
1-46 -4 1 Ya w (23) 
1 : 
si 1. -4 6 -4 1 Yi = w(zi) |, 
Ll. 4 q —£ 1 YN-3 w (Ty-3) 
1 dB A UYN-2 w(2y-2) 
1 4 7 UN—1 w(ty-1) 
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2.4 Métodos de una etapa. Ejemplos 


de discretos continuos 
h= f= 0.5 | (£.=.0.25 

0.00 0.0000 | 0.0000 0.0000 0.0000 
0.25 2.1178 2.1013 
0.50 4.2862 4.1873 4.1546 
0.75 6.1611 6.1126 
1.00 8.9795 | 8.1861 7.9938 7.9302 
1.25 9.6429 9.5692 
1.50 11.3445 | 11.0700 10.9793 
1.75 12.2413 12.1390 
2.00 | 14.8767 | 13.4689 | 13.1289 13.0167 
2.25 13.7120 13.5916 
2.50 14.3603 | 13.9773 13.8509 
2.15 13.9204 13.7904 
3.00 | 15.6096 | 13.9450 | 13.5468 13.4156 
3.25 12.8725 12.7430 
3.90 12.3053 | 11.9254 11.8001 
3.19 10.7441 10.6253 
4.00 | 11.0959 | 9.7111 9.3770 9.2667 
4.25 7.8811 7.7811 
4.50 6.5886 6.3214 6.2332 
4.75 4.7698 4.6946 
5.00 4.2534 | 3.4889 3.3038 3.2427 
5.25 2.0061 1.9597 
9.90 1.0569 0.9629 0.9319 
9.19 0.2638 0.2483 
6.00 0.0000 | 0.0000 0.0000 0.0000 


Tabla 2.7: Resultados obtenidos, By, para distintos valores de h y solución exacta. 


que podemos utilizar para resolver con una retícula de paso h = 0.5 (figura 2.15(c)) así como 
para una red de paso h = 0.25 con los nodos numerados como en la figura 2.15(d). 

Las gráficas asociadas al conjunto de valores hallados para cada paso de red pueden ser 
comparadas entre sí y con los valores exactos en la figura 2.16. Los resultados para los 
distintos valores de h pueden compararse en la tabla 2.7. 

Eligiendo los resultados obtenidos para la proyección sobre la red Ry-0.25, pueden calcularse 
los momentos flectores con las fórmulas: 


— Yi+1 — 2Yi + Yi-1 
mi > 
0.25% 


cons=0,1.2. 0 2% 


Se completan las funciones de deformación y esfuerzo obteniendo los giros de las secciones 
con auxilio de las fórmulas: 


ia AL, con2¿=0,1,2,...,24, 
0.5 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


16 16 16 
14 14 14 
12 12 12 
10 10 10 
8 8 
6 6 
4 4 
2 exacta 2 exacta exacta 
—6— h=1 —6— h=0.5 —0—h=0.25 
] : ñ : ] 0 : j j á j e : j j ] ñ 
0 1 2 3 4 5 6 0 1 2 3 4 5 6. 0 1 2 3 4 5 6 


Figura 2.16: Solución analítica junto con la solución numérica para h= 1, h=0.5 y h = 0.25. 


y finalmente, para obtener los esfuerzos de corte aplicaremos: 


e Via = Liar + ¿ia == Via 


i 2.0253 , con2=1,2,...,23. 


La ley de cargas, los esfuerzos de corte, los momentos, los giros de las secciones y las 
flechas halladas para la aproximación h = 0.25 pueden verse representadas por las gráficas 
de la figura 2.17. Los descensos y giros deben ser divididos por el coeficiente B. 


4 T T J T T 
0 1 Jl 1 1 L 
0 1 2 3 4 5 6 
10 T T J T T 
E oP AA J 
-10 1 j i i j 
0 1 2 3 4 5 6 
10 T T T T T 
= 0 | 
-10 i j i i j 
0 1 2 3 4 5 6 
10 T T T T T 
s 0 
=10 1 l 1 1 L 
0 1 2 3 4 5 6 
20 T T J T T 
a 101 ] 
0 i j i i | 
0 1 2 3 4 5 6 


Figura 2.17: De arriba a abajo: ley de cargas By?”, cortante By””, flectores By”, giros By” y flechas By. 
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2.4 Métodos de una etapa. Ejemplos 


2.4.6. Viga doblemente empotrada 


Viga de 8 m doblemente empotrada sometida a una ley de cargas, simétrica respecto al 
eje vertical que pasa por su punto medio, de envolvente parabólica definida por la función 
w(x) =—a? +8x (figura 2.18(a)). 

Podemos establecer una red de paso h = 1 m. Las condiciones de simetría nos permiten 
numerar los nodos como en la figura 2.18(b) y reducir el sistema de ecuaciones a la mitad. 


(a) 
16t 


6 x=4 
x=0 lI=8m x=8 


(b) 


15 16 


5 15.75 16 


13.75 1 
2 


(d) 


2 8 1543 15:75 193 16 


Figura 2.18 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


x discretos continuos 
h=1 h=0.55|h=0.25 

0.00 0.0000 0.0000 0.0000 0.0000 
0.25 2.1333 2.0223 
0.50 8.9332 7.8673 7.6465 
0.75 16.5435 16.2153 
1.00 34.1250 | 28.8203 | 27.5181 27.0861 
1.25 40.1685 39.6373 
1.50 99.1832 | 53.8995 93.2746 
1.75 68.1490 67.4366 
2.00 94.5000 | 84.7813  82.3926 81.6000 
2.25 96.1490 95.2839 
2.50 111.7832 — 108.9835 108.0545 
2.15 120.5122 119.5281 
3.00 || 146.1250 | 133.5078 | 130.4048 129.3750 
3.20 138.3872 137.3216 
3.90 147.5332 | 144.2443 143.1527 
3.19 147.8208 146.7137 
4.00 | 166.0000 | 152.3750 | 149.0234 147.9111 


Tabla 2.8: Resultados obtenidos, By, 


Se puede escribir el sistema: 


para distintos valores de h y solución exacta. 


y-1— Yo +6y1 —4yat+ ya = 7, con Y-1 = Ya, Ya =Ú 
Yo — 4y1 + 6ya — Aya + Ya = 12, con Yo =0 
yi — Aya + Óya — Aya + yg = 15, 
Ya —Ays + Oy d+ y = 16, 
que matricialmente es: 
7 —d4 1 0 Yi i 
-4 6 -4 1 y | | 12 
1. -4 7 -4 y | |15 
0 2 -8 6 Ya 16 
y cuya solución es: 
Yo 0.0000 
Yi 34.1250 
ya | = 94.5000 
Y3 146.1250 
Ya 166.0000 
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2.4 Métodos de una etapa. Ejemplos 


exacta 20 
—6— h=1 


exacta 4 20 
—8— h=0.5 


exacta 
—8— h=0.25 


0 2 4 6 8 


0 2 4 6 8 0 2 4 6 
Figura 2.19: Solución analítica junto con la solución numérica para h= 1, h=0.5 y h = 0.25. 


Para un valor de h genérico, el sistema se puede escribir como: 


A Yi w (21) 
4 6 -4 1 
1. -4 6 -4 1 e a 
Y3 w (23) 
1 o en ade E 2d 
7 1-4 6 -4 1 gy |=| vwíx) |, 
1-4 6 -4 1 UN-2 w(ty-3) 
1.A<A Y == YN-—1 w(2y-2) 
2 38 6 YN w(en-1) 
donde en este caso, N = qa 


Utilizando esta última expresión dada en forma genérica, resolvemos para la configuración 
de la figura 2.18(c), habiendo tenido en cuenta una vez más las condiciones de simetría y 
para una red de paso h = 0.5. Asimismo, resolvemos también para la configuración de la 
figura 2.18(d), esta vez para una red de paso h = 0.25. 

Los resultados para los distintos valores de h pueden ser comparados entre si y con los 
valores exactos en la tabla 2.8 y pueden también visualizarse las aproximaciones obtenidas 
en las tres redes sobre la gráfica de la figura 2.19 (donde además se ha dibujado la solución 
aproximada sobre todo el intervalo aprovechado la simetría que presenta la solución). 

Eligiendo los resultados obtenidos para la proyección sobre la red Ry=0.25, pueden ser cal- 
culados ahora los momentos flectores en sus puntos a través de los cálculos complementarios: 

mi HE E o O 
donde debe entenderse que Y¡7 = Y15. 

Podemos completar las funciones de deformación y esfuerzo obteniendo ahora los giros de 
las secciones con: 


yl = Yi41 7 Yi-1 

] 0.5 ! 

nuevamente con Y17 = Y15, y para obtener, finalmente, los esfuerzos de corte aplicaremos: 
m Viso — 2Yiya + 241 — Yi=2 


mM — o E 
Yi 20.253 y 


con ¿=0,1,2,...,16, 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


Figura 2.20: De arriba a abajo: ley de cargas By”, cortante By””, flectores By”, giros By” y flechas By. 


donde ahora Yi7 = Y15 Y Y1s = Y1a: 

La ley de cargas, los esfuerzos de corte, los momentos, los giros de las secciones y las flechas 
halladas sobre la mitad izquierda de la viga para la aproximación h = 0.25 pueden verse 
sobre las gráficas de las figura 2.20. Para poder completar las gráficas en todo el intervalo se 
ha tenido en cuenta que las derivadas de orden par son funciones pares con respecto el eje 
vertical de simetría y las impares son funciones impares. 


2.5. Los métodos de dos etapas en ausencia de cargas 
y momentos aislados. Ejemplos 


Como alternativa a los problemas resueltos anteriormente y como introducción a los méto- 
dos que utilizaremos sistemáticamente en el $ 3 podemos plantear ahora la integración de la 
ecuación de cuarto orden (2.2) en dos etapas, cada una de las cuales aborda una ecuación 
de segundo orden. En la primera etapa se obtendrían los momentos flectores: 


m” (2) =w (2), 
en tanto en la segunda obtendríamos las flechas que satisfacen la ecuación: 
By (a) =milzx). 


La integración aproximada de ambas ecuaciones simplifica en cierta medida el algoritmo 
aunque multiplica, en general, el número de ecuaciones ya que ambas etapas conducirán ge- 
neralmente a un solo sistema de ecuaciones, es decir: los momentos sobre los nodos quedarán 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


tipo de viga extremo izquierdo extremo derecho 


Y. 


ménsula empotrada en su 
extremo izquierdo 


a 


ménsula empotrada en su 
extremo derecho 


2 A 


viga doblemente apoyada 


Ha 


viga empotrada-apoyada 


A 


viga apoyada-empotrada 


J | 


viga doblemente empotrada 


Tabla 2.9: Leyes de cargas segmentariamente continuas sin acciones aisladas en los extremos. 


normalmente ligados a los valores de las flechas. Los dos sistemas de ecuaciones algebraicas 
podrán ser resolubles sucesivamente, en dos etapas reales, solamente en el caso de ménsulas 
o vigas doblemente apoyadas (casos isostáticos), en tanto que ambos deberán ser resueltos 
simultáneamente en el resto (casos hiperestáticos). 

Las condiciones de extremo en el caso en que las leyes de cargas no incluyan cargas 
concentradas ni momentos aislados deben ser reinterpretadas en términos de momentos y de 
flechas lo que reduce los puntos ficticios al primer nivel, es decir: a los puntos a distancia 
h, del extremo. Estas condiciones vendrán ahora dadas en términos de las funciones (flechas 
y/o momentos) o de sus derivadas primeras según la tabla alternativa 2.9. 

Las condiciones de extremo tipo deben ser asimismo replanteadas a la hora de asignar 
valores a los puntos “ficticios”. Tomando el extremo izquierdo xy y el nodo interior 1, como 
referencia, se muestran en la tabla 2.10 las condiciones de extremo y los valores que se 
deben asignar con este procedimiento al punto ficticio _, según las condiciones del vínculo 
extremo respecto al valor del punto interior próximo al punto frontera 2y. Estas condiciones 
se completan en la tabla con yo = 0, descenso nulo en apoyos rígidos y empotramientos 
perfectos y my = 0 en apoyos. Estos resultados son válidos en el extremo derecho sin más 
que cambiar la notación de los subíndices. 

Para ilustrar el método en "dos etapas" se retoman tres de los ejercicios anteriores sobre 
una red de paso h = 0.25 m y se obtienen directamente los momentos flectores y las flechas 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


vínculos condiciones de extremo 
1% etapa 2% etapa 
empotramiento perfecto Yo =0 Yo =0 


% y 
Ó 1 


FP—_—_—. 0 yy 0 = 
2h Ya = Yi 
apoyo rígido my =0 Yo =0 
O 
extremo libre mo =0 mo =0 
4-4 + 
e] o 1 m3 = M1 0 
o m_¡ =m 
2h 1 1 


Tabla 2.10: Solicitaciones segmentariamente continuas sin acciones aisladas en los extremos. 


de cada pieza. 


2.5.1. Ménsula empotrada en su extremo derecho 


Corresponde al ejemplo descrito en $ 2.4.2; se trata de la viga de 4 m en libre-empotrada 
cuya ley de cargas venía definida por la función: 
e 4x1, si0<x=<, 
212 “sl<z<4, 
según la figura 2.6(a). 


Sobre el nuevo esquema de la figura 2.21 que incluye un solo punto exterior a cada uno 
de los lados del intervalo que determina la longitud de la pieza, plantearemos el sistema de 


Figura 2.21 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


ecuaciones: 
12 etapa de etapa 
02 (M1 —2mp+ mx) = 0, 2 ly 2Y1 +90) = m, 
035 (Ma — 21 + mg) =, L, E (yz — 2Ya +Y1) = UE 
0257 (M3 — 2m93 + m1) = 2, E 7 (Ya 2Y3+Y2)  = Ma, 
ES (my — 2m3 + ma) == me 7 (Ys — 2Ya + Ys) = Ma, 
537 (ms — 2m4 + m3) = A, me > (Yo — 245 + Ya) = ¡Ms 
537 (mó — 2m5 + m4) = 0, ne (yr — 2Y6 + Ys) = M6, 
13 (m7 — 2m6 + m5) = 6, O. 553 (Ys — 2y7 + Y6) = M7, 
025 (mg — 2m7 + mg) = 03 > (Ya — 2Yg + yr) = Mg, 
Ae (mo — 2m3 + mr) = 0, E 2 (Y1o — 2Y9 + ys) —= Mg, 
om (mio — 2mM9 + mg) = 2, E 2 (Yi == 2410 ES Yo) = Maio, 
025 (Mi — 2M10 + Mg) = 2, Me > (Y —2Y1 +Y10) = Mar, 
535 (mi — 211 + M10) = 2, 03% 2 (Ys — 2Y12 + Yi) = Mi», 
025 (mi — 2m12 + M1) = 2, > (Yia — 2Y13 + Yi2) = Maz, 
035 (Mis — 293 + Mia) = 2 lus 2Ya+Ys) = Ma, 
035 (Mis — 2M14 + Mig) = 2, E 2 (Yo — 2Y15 + Ya) = M5, 
02 (Mis — 22M + Mia) = 2 lu Y +Ys) = Ms, 
con las condiciones: con las condiciones: 
mo =0, Yy16 =0, 
m-1= my, Yi? —= Y15- 


Este sistema contiene dos subsistemas de Cramer, es decir: resoluble en dos etapas suce- 
sivamente, en la primera de las cuales podemos hallar los momentos y posteriormente un 
segundo sistema dependiente de los valores hallados en el primero. 

Las soluciones de momentos y flechas del sistema completo son: 


mo 0.0000 Yo 98.7031 
ma 0.0000 YN 90.7500 
ma 0.0625 Ya 82.7969 
ma 0.2500 Ya 74.8477 
ma 0.6250 Ya 66.9141 
ms 1.2500 Ys 59.0195 
me 2.1875 Yo 51.2031 
ma 3.5000 yr 43.5234 
ms |=| 5.2500 |; ys |=| 36.0625 
mao 7.3125 Yo 28.9297 
m0 9.5000 Y1o 22.2539 
ma 11.8125 Y11 16.1719 
m9 14.2500 Y13 10.8281 
m3 16.8125 Ys 6.3750 
Ma 19.5000 Ya 2.9727 
mas 22.3125 Ys 0.7891 
mi 25.2500 Y16 0.0000 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


No se debe olvidar que los valores obtenidos para la función y deben de ser divididos por 
el factor de rigidez B con el fin de que proporcionen el valor de los descensos en su verdadera 
magnitud. Pueden compararse estos resultados con los obtenidos para momentos y flechas 
en el ejemplo descrito en $ 2.4.2 comprobando que son resultados idénticos. 


Para facilitar la resolución del problema, podemos escribir en forma matricial para un 
valor genérico de h: 


donde 
2 ma w (29) 
DH ma w (21) 
7 il — 1 m3 w (22) 
1 —2 1 MN-1 (00) (Ty-2) 
1 —=2 1 MN (00) (2y-1) 
y 
1-2 1 Yo 
1 2 1 Yi 
1 ; Ya 
A) = 72 L == A E Ya 
1-2 1 : 
1 -2 YN-2 
2 YN-1 


2.5.2. Viga doblemente apoyada 


Corresponde al ejemplo descrito en $ 2.4.3; se trata de la viga de 6 m. con dos apoyos 
rígidos en sus extremos solicitada por una ley de cargas que puede venir representada por la 
función: 


—2+4, si0<zx<2, 


wta=x4 1 si2<r<L 
=x+ Y, si4<x<6. 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


Sobre el nuevo esquema de la figura 


123450672 


el sistema de ecuaciones: 


1% etapa 

037 (m2 — 2m1 + mg) 
0.257 (m3 — 2m9 + m1) 
0.257 (my — 2m3 + mo) 
0.35 (ms — 2m4 + m3) 
037 (m6 — 2m5 + ma) 
037 (m7 — 2m6 + ms) 
037 (mz — 2m7 + mg) 
037 (my — 2m3 + my) 
037 (mio — 2m9 + mg) 
037 (mi1 — 2m10 + Mg) 
0.257 (mia — 2m11 + Mio) 
033 (mi3 — 2m12 + m1) 
0.957 (Mi, — 2m13 + M9) 
0.35 (mis — 2M14 + Mi) 
0.357 (mis — 2m15 + M4) 
033 (mi7 — 2m16 + mi5) 
037 (mis — 2m17 + M6) 
037 (mi9 — 218 + M7) 
037 (Mon — 2m19 + Mig) 
037 (ma, — 2m90 + Mig) 
0.957 (Mo, — 221 + Map) 
037 (ma — 2m92 + Ma1) 
0.252 (Mo, — 293 + Mao) 


con las condiciones: 
mo = 0, 


11124 


=0, 


175 15 


125 1 


9 10 11 12 13 14 15 16 17 12 19 20 11 22 233 


Figura 2.22 


2d 


2.22, sin puntos de red exteriores, plantearemos ahora 


eN etapa 
0357 2 (Ya — 2y1 + Yo) 
0.357 3 (Y — 2Y2 + Ya) 
025 (Ya — 2Y3 + Ya) 
dE 5 (Ys — 2Ya + Ya) 
e > (Yo — 2Y5 + Ya) 
025 (Y7 — 2Y6 + Ys) 
de 3 (Ys — 247 + Ye) 
025 (Yo — 2Ys + yr) 
025 3 (Y1o — 2Yg + Ys) 
sde 7 (Y11 — 2Y10 + Yo) 
Y 2 (Y12 = AU + y10) 
ale 7 (Y1g — 2Y12 + Y11) 
575 7 (Y1a — 2Y13 + Y12) 
Sos 3 (Y1s — 2Y14 + Y13) 
E 7 (Y1ó — 2Y15 + Y1a) 
e 7 (Y17 — 2Y16 + Ys) 
Y 3 (Yis — 2Y17 + Y1s) 
538 7 (Y19 — 2Y1g + Y1r) 
025 7 (Yoo — 2Y19 + Ys) 
Y 3 (Yo — 2Y20 + Y1o) 
Y 7 (Ya — 2Y21 + Yao) 
0 7 (Yoz — 2Ya22 + Ya1) 
$ 52 (Yaa — 2Yo3 + Yoo) 


con las condiciones: 


Yo = 0, 
ya, = 0. 


La primera etapa es un sistema de Cramer resoluble independientemente. Ambas etapas 
constituyen por lo tanto sistemas resolubles sucesivamente. Podemos hallar en primer lu- 


gar los momentos y posteriormente las flechas o resolver el sistema de momentos y flechas 


simultáneamente. 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


De forma matricial, y para un valor de » genérico, los sistemas los podemos escribir como: 
AM =F; AY = M 


donde en este caso: 


=2 1 
1 -2 1 
1 1 -—2 1 
A= h? ñ Eo ) 
1-2 1 
1 -—2 

mi Yi w (21) 

Ma Ya w (22) 

mel E de Bel le pel? (23) 

MN-2 YN-2 w (IN-2) 
MN-1 YN-1 w (un-1) 
Para h = 0.25, las soluciones son: 

mo 0.0000 Yo 0.0000 
mi —1.5990 Yi 4.0677 
ma —2.9635 Yo 8.0354 
ma —4.1094 Y3 11.8179 
Ma —5.0521 Ya 15.3435 
ms —5.8073 Y5 18.5534 
Mé —6.3906 Y6 21.4004 
m; —6.8177 Y7 23.8479 
Ma —7.1042 Us 25.8694 
My —7.2969 Y9 27.4468 
mia |= —7.5000 |; ya | =| 29.4141 
M13 —7,4427 Y13 29.1341 
Mia —7.3229 Ya 28.3889 
m5 —7.1406 Y15 27.1860 
Mi6 —6.8958 Yl6 25.5369 
m7 —6.5260 Y1i7 23.4568 
Mi8 —5.9844 Ys 20.9687 
m9 —5.2865 Y1i9 18.1067 
Mo 2.4115 Y22 7.7556 
Mag —1.2448 Y23 3.9167 
11194 0.0000 Y24 0.0000 


que coinciden con las obtenidas en el ejemplo descrito en $ 2.4.3. 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


2.5.3. Viga apoyada-empotrada 


Corresponde al ejemplo de $ 2.4.5; se trata de una viga de 6 m con un apoyo rígido en 
su extremo izquierdo y un empotramiento perfecto en el derecho solicitada por una ley de 
cargas que puede venir representada por la función: 


—zx +6, 


si0<zx<3, 
si3<x=<6. 


Sobre el nuevo esquema de la figura 2.23, con solamente un punto ficticio al otro lado del 
empotramiento, plantearemos el sistema de ecuaciones: 


035 (ma 2m1 + ma) 
035 (m3 — 2m3 + m1) 
0.357 (my — 2m3 + mo) 
0.357 (m; — 2m4 + m3) 
037 (m6 — 2m5 + ma) 
033 (m7 — 2m6 + ms) 
037 (mg — 2m7 + mg) 
035 (my — 2m3 + m+) 
037 (mio — 2m9 + mg) 
0.35 (mi1 — 2m10 + Mg) 
0.357 (mia — 2m11 + Mio) 
037 (mi3 — 212 + m1) 
0.35? (Mia — 2m13 + M9) 
0.35? (mis — 214 + M13) 
0.35" (Mis — 2M15 + Maa) 
033 (m7 — 2m16 + mi5) 
037 (mis — 2m17 + M6) 
02357 (Mi9 — 218 + M7) 
0357 (M2 — 219 + Mag) 
037 (ma — 2M20 + Mag) 
0,37 (m2 — 221 + Map) 
037 (ma — 292 + Ma1) 
027 (Mas — 2M233 + Maz) 


con la condición: 
mo = 0, 


Y2as — 2Y23 + Y22 
Yos — 2Yo4 + Yoz 


027 (Ya — 2Y1 + Yo) 
pe 3 (Y3 — 2y2 + Y1) 
e (Ya e 2Y3 + Ya) 
so (Ys — 2ya + Y3) 
035 (Ys — 245 + ya) 
oa (Y7 — 2Y6 + Ys) 
025 (Ys — 2Y7 + Y6) 
0357 (Yo — 2yg + yz) 
027 (Y1o — 2Yo + Ys) 
0. 538 (Yy11 — 2Y10 + Yo) 
0257 (Y12 — 2Y11 + Yo) 
02 (Y13 — 2Y12 + Y11) 
0257 (Ya — 213 + Yr2) 
o. 553 (Yis — 2Y14 + Yiz) 
077 (Y16 — 2Y15 + Ya) 
0. 53 (Yir — 2Y16 + Yis) 
O. 55 (Y1is — 2Y17 + Yis) 
a (Yi — 2Y1g + Ya7) 
538 (Yao — 2Y19 + Ys) 
a (Y2a1 — 2Y20 + Y1o) 
0357 (Y22 — 2Y21 + Yoo) 
Ñ y (Yo — 2Y292 + Ya1) 
1 
( ) 
7 l ) 


Yo =U, 
Y2s = 0, 
Yas — Y23 


Las dos etapas de este sistema deben ser resueltas simultáneamente. El conjunto de ambas 
etapas constituye aquí un solo sistema de Cramer. 
Escrito matricialmente para un A genérico tendremos: 


AX =G, 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


Figura 2.23 


donde ahora la matriz A se puede escribir como una matriz por bloques, así como los vectores 
X y G, de la siguiente manera: 


siendo /y la matriz identidad de dimensión Nx N, 0 denota una matriz (o un vector) con to- 
dos los elementos nulos, mientras que las matrices Aj €  Miw-uyxvw (R) y 
Ao € Mnx(w 1) ( R) son: 


=2 1 
E ll =Y 1 
1 2 1 1. -2 1 
1 -—2 1 
A o ; An => na 
1 -2 1 , ES a 
1 -2 1 ) 
y los vectores M, Y y F": 
ma Yi wW (21) 
ma Ya 
m3 Y3 E (2) 
M= ! Y = : F = w (13) 
MoNn-1 UN-2 
wW(Ty- 
Mon UN-1 ( Al 1) 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


Las solución de este sistema para h = 0.25 es: 


Mo 0.0000 YO 0.0000 
mi 0.7733 Yi 2.1178 
ma —1.5309 Ya 4.1873 
ma 335% Ys 6.1611 
Ma 2.9368 Ya 7.9938 
ms 3.9938 Y; 9.6429 
m6 —4.0927 Y6 11.0700 
mx —4.5378 yr 12.2413 
Mg —4.8739 Y8 13.1289 
Mo 5.0843 Yo 13.7120 
Mio —5.1544 Y10 13.9773 
11111 —5.0683 Y11 13.9204 
mi |=| -4.8103 |; ya |=| 13.5468 |, 
M13 —4.3648 Y13 12.8725 
mi15 2.9738 Y15 10.7441 
Mi6 2.0596 Y16 9.3770 
mMai7 —1.0203 Y17 7.8811 
M8 0.1283 Yis 6.3214 
Mig 1.3707 Y19 4.7698 
1130 2.6912 Y20 3.3038 
1Ma1 4.0742 Ya1 2.0061 
Maa 9.5040 Y22 0.9629 
Mas 6.9651 Y23 0.2638 
Maa 8.4419 Y24 0.0000 


que es la obtenida en el ejemplo de $ 2.4.5. 


2.5.4. Viga doblemente empotrada 


Corresponde al ejemplo de $ 2.4.6; se trata de la viga doblemente empotrada de 8 m 
con condiciones de simetría; está solicitada por la ley de cargas de envolvente parabólica 
w(x) =—a2? + 8x. 


Sobre el nuevo esquema de la figura 2.24, con solamente un punto ficticio en el empo- 
tramiento izquierdo, habiendo tenido en cuenta las condiciones de simetría en la numeración 
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2 Solicitaciones segmentariamente continuas con saltos finitos 


14.43 15 15.4315.7515.93 16 


Figura 2.24 


de los nodos, el sistema de ecuaciones será: 


> (ma — 2m1 +mo) = 19375  ia(Y-2M0+Yy1) = Mm, 
0357 (m3-2m2+m) = 3.75, 57 (yo — 2y1 + Yo) ==" a 
035 (ma —2m3+ma) = 5.4375, 57% > (ys — Za + y) => as 
075 (Ms — 2M4 + M3) = A (Ya — 243 + Ya) = mM, 
a (mé —2m5+m4) = 8.4375, 03 > (Ys — Zya + ys) = ma, 
032 (M7 — 2mM6 + M5) = 9.75, aL 2 (Yo — 2Y5 + Ya) => MES 
(ma —2m7 + me) = 109375, oi(yr—2Y0+Y)  = Mo 
> mo 12, 0. 27 (Ys — 2y7 + Ys) = m7, 
035 (Mio — 2m9 + Mg) = 12.9375, — 5a5x (Yo — 2Ys + yr) = Mg, 
0. 027 (m1 —2mMi9+ mg) = 13.75, a) 27 (Yo —2Y9 + Ya) = Mo, 
035 (mi — 211 + M0) = 14.4375, 0. 7 (Ya — 2Y10 + Yo) = Mio, 
037 (ma3 — 212 + Mu) = 15, 0. 7 (Y12 — 2Y11 + Yio) = Maz, 
035 (Mia — 213 + Ma2) = 154375, els Ye+Yu) = Maz, 
537 (Mis — 214 +mM13) = 17.75, 0. 27 (Yia — 2Y13 + Y12) = Maz, 
557 (mi — 2m15 + ma) = 15.9375, (2 +93) = ma 
03 (Mas — 2M6 + Mas) = 16, 027 (Y16 — 2Y15 + Yia) = Ma, 
027 (Y1s — 2Y16 + Yis) = Ma, 
con las condiciones: 
Yo =0, 
Y-1= Y1- 


Ambas etapas deben ser resueltas simultáneamente. Escrito de forma matricial para un 
valor de h genérico, el sistema se puede escribir como 


AX =G, 


donde la matriz A se puede escribir como una matriz por bloques, así como los vectores X 
y G, de la siguiente manera: 


Ar 0 M F 
A= A X=Y|—|]: G=|- 
Gasta): *= (7): e- (5) 


siendo [y +1 la matriz identidad de dimensión (N + 1)x (WN + 1) (en este caso N = 5% debido 
a la simetría), O denota una matriz (o un vector) con todos los elementos nulos, mientras 
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2.5 Los métodos de dos etapas en ausencia de cargas y momentos aislados. Ejemplos 


que las matrices A11 € Myx(w+1) y A22 € Miv+nxn son: 


mo Yi w (11) 
mi Ya w (23) 
We ee Yo Y3 FP= w (23) 
MN-1 YN-1 w(2y-1) 
MN YN w(2y) 
La solución para h = 0.25 es: 
Mo 68.2666 Yo 0.0000 
mi 57.6104 Yi 2.1339 
ma 47.0752 Ya 7.8673 
m3 36.7744 Y3 16.5435 
Ma 26.8135 Ya 27.5181 
m5 17.2900 Ys 40.1685 
M6 8.2939 Y6 53.8995 
m7 —0.0928 Y7 68.1490 
Ma = —7.:7959 : Ys = 82.3926 , 
Mg —14.7490 Yo 96.1490 
Mi0 —20.8936 Yio 108.9835 
M1 —26.1787 Yi 1200122 
mM12 —30.5615 Y12 130.4048 
M13 —34.0068 Y13 138.3872 
Mia —36.4873 Y1a 144.2443 
Mi5 —37.9834 Yi5 147.8208 
Mi6 —38.4834 Y16 149.0234 


que es la misma que la obtenida en el ejemplo de $ 2.4.6. 
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